Robustness and diffusion of pointer states 
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Classical properties of an open quantum system emerge through its interaction with other degrees 
of freedom (decoherence). We treat the case where this interaction produces a Markovian master 
equation for the system. We derive the corresponding distinguished local basis (pointer basis) by 
three methods. The first demands that the pointer states mimic as close as possible the local non- 
unitary evolution. The second demands that the local entropy production be minimal (predictability 
sieve). The third imposes robustness on the inherent quantum and emerging classical uncertainties. 
All three methods lead to localized Gaussian pointer states, their formation and diffusion being 
governed by well-defined quantum Langevin equations. 



The programme of decoherence has been very successful in explaining the classical appearance of macroscopic or 
mesoscopic quantum systems, both theoretically and experimentally [jp. The interaction of a quantum system with 
its environment leads in these cases to a derealization of phase relations in the full configuration space of system plus 
environment, preventing them to be observed locally, i.e. at the system itself. Thereby the environment distinguishes 
a certain preferred basis for the system, which can be used to describe the apparent classical behavior (pointer basis 
Q). An important question is how the pointer basis can be determined and whether it is uniquely fixed. It would 
then be possible to give a unique decomposition of the reduced density matrix into an apparent ensemble of wave 
functions. 

No unique rules have so far been adopted to calculate the pointer states in the general case. In [| the suggestion 
was made that the pointer basis (there called collection of "memory states" ) is characterized by its robustness (there 
called "dynamical stability" ) . A first quantitative measure to investigate the dynamical stability is the rate of de- 
separation introduced in [Q — it measures how fast a quantum system becomes entangled with environmental degrees 
of freedom. In a model consisting of harmonic oscillators, it was shown that coherent states are the most stable states 
and therefore can be considered as pointer states |l|, Q . A different measure for robustness is the "predictability sieve" 
put forward in The pointer basis is there distinguished by the property of having the least production rate for 
local entropy during the coupling to the environment. In the case of harmonic oscillators, this again leads to the 
coherent states as the pointer basis At least for such simple systems, rate of de-separation and predictability sieve 
are roughly equivalent measures Q. 

On the other hand, the theory and formalism of quantum state diffusion (QSD) were put forward to attribute 
random wave functions for local systems, which satisfy an appropriate Langevin equation 0. These wave functions 
are known to be related to possible continuous measurements (8) as well as to decoherent histories Q of the given 
local system. But even if we take them as wave functions of mere formal meaning (since subsystems do not, in 
general, possess their own pure states), the question arises whether there is any connection between these states and 
the pointer states, in cases where the local system exhibits classical properties. We shall show that there is in fact 
such a connection - pointer basis and QSD basis are substantially the same. For this aim, we shall also present below 
a new, alternative, derivation of QSD. 

In the following we shall consider the dynamics of the reduced density matrix, p(t), of a system interacting with a 
certain decohering environment. Ideal pointer states (described by a fixed set of projectors P n ) are characterized by 
the fact that p(t) can be decomposed as 

P(t) -> f»P» ' t > tD ' W 
n 

for a generic initial state p{0), where to is the decoherence time. The weights {f n } correspond to a normalized 
probability distribution. The pointer states {P n ] H = 1, 2, . . .} form in this case an orthogonal system. For macroscopic 
systems, tr> is extremely short Jl], px[ . More generally, one would expect 

p(t) -> J f(T)P(T)dT , t » t D , (2) 
where f(T) is a probability distribution over the pointer states P(T) which project now on an overcomplete set of 
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pure states (the above-mentioned coherent states provide an example for this). The pointer states in ([j]) or (||) result 
after an explicit interaction with the environment is taken into account |Q . 

In many cases, the effects of decoherence can be described by the following Markovian master equation H, 

^L^Cp{t) = -^-\p2,p(t)]-j[x,[x,p(t)}} , (3) 

where D describes the strength of the interaction with the environment. Such an equation arises, for example, in 
situations where environmental degrees of freedom scatter off a macroscopic object and localize it by carrying away 
quantum correlations with the object jl], Applying the concept of predictability sieve would mean to minimize the 
local production of "linear entropy" S(t) = 1 — tr p 2 (t). This does not give a unique answer, since the result depends 
explicitly on t. In the oscillator case, one has therefore calculated the time-integrated entropy production [jjj §]. If 
one considers the initial entropy production rate 5(0), assuming the initial state p(0) of the subsystem to be a pure 
state P(r), one finds from (|) that 5(0) = D ((x 2 ) - (x) 2 ) = Da 2 , where the expectation value refers to the initial 
state. This rate would be minimized if the pointer wave functions were delta functions. However, their spread a 
increases dynamically due to the kinetic term in (Q) . Therefore, very narrow wave functions do not produce minimum 
entropy on a finite time scale and thus cannot correspond to pointer states. The (coherent) unitary spreading and 
the (incoherent) non-unitary localizing terms of the master equation (|3|) are competing with each other. For a wave 
function of characteristic width er, the above two effects are approximately balanced for the "equilibrium width" 

a ~ (Cm)" 1/4 . (4) 

It is then reasonable to conjecture that, in the spirit of predictability sieve and of dynamical robustness, <jq will be 
the characteristic width of the pointer states. This is in fact what we shall show by using three different methods, all 
invoking a principle of robustness. 

The first method goes as follows. Let us allow for the pointer state P(T) a certain natural time dependence such 
that it may initially evolve as close as possible along the true state p satisfying the master equation (^), and then 
reach a stationary state. We introduce the "speed" v describing the departure of the states P(T) from p in the 
Hilbcrt-Schmidt norm: 



v 2 = tr 



d 



dt P(T) - CP(T) 



(5) 

The smaller v, the greater is the robustness of the pointer states P{T). Hence one defines the optimum drift of the 



pure pointer state P(T) = ipip' by minimizing v [p2[. This is given by the nonlinear equation [112, 13 



ii = (A/^ f )V> - (c^)tp . (6) 

This result is valid for all kinds of Markovian subdynamics. In our special case ([|), it yields the following nonlinear 
wave equation: 

i> = -^fi> - f [{£- <x>) 2 - a 2 ] i> . (7) 

As shown in pL this equation has a stationary solution which is unique up to Galilean transformations. The wave 
function of the fiducial stationary state is the complex Gaussian wave packet 

Ms) = («fl/27r) 1/4 exp(-ax 2 /4) , (8) 

with parameter 

a = a R + iai = (1 - i)V2Dm . (9) 

The principle of "Hilbert-Schmidt robustness" has thus singled out unique Gaussian pointer states as the robust pure 
states closest to the true non- unitary local dynamics. The exact width a = \j Jur confirms the heuristic estimate 
(Q). Accordingly, we restrict our further discussion to pointer states P(T) with Gaussian wave functions and make 
the ansatz 

ip r (x) = (a R /2ir) 1/ ' l exp(-a(x-x) 2 /4 + ip(x-x)) , (10) 



3 



where T = (x,p) T has been understood. For later purposes, we calculate the correlation matrix C 



(xp + h.c.)/2 



(xp + h.c.)/2 
p 2 



l^o) = — 



1 - ai /2 
-aj/2 |a| 2 /4 



(11) 



of the quantum uncertainties of the canonical observables in the pointer state itself, where "0o denotes the fiducial state 
(||). It was shown in [ |To| that the states diagonalizing p exactly are the harmonic oscillator eigenfunctions which are 
very broad, while narrow eigenfunctions are apparently obtained only for discrete systems. In contrast to these, the 
above pointer states are well localized. For example, in the situation of a small dust particle (to = 1CP 14 g) scattered by 
air molecules one has D ~ 10 32 cm _2 s _1 and therefore do ~ (Z?m) -1 / 4 w 10~ n cm and trj ~ \Jm/D « 1CP 10 s. 

We now come to the second method. As a preparation, we shall discuss the reduced dynamics of the local system 
in the basis given by (|Io|). We allow temporarily the parameter a to take an arbitrary complex value, and then derive 
again a distinguished value. If one allows a "natural" time dependence for the probability distribution f(T; t) of the 
pointer, the asymptotic condition (0) can be turned into an exact identity: 



p(t) = J f(T;t)P(T)dT, t>t D , 



(12) 



where dr = dxdp/27r. This important fact will be proven elsewhere |lq ]. It generalizes the corresponding statement 
made in |l4| for the specific value ([)]) of a as well as the asymptotic statement proved in |l^, [TtJ . 
From (|) and (|l|) one can derive an evolution equation for f(T;t): 



df(T;t) p 



d s f(T; t) + - [D pp d 2 pp + IK . Ut + 2D px d 2 x ] f(T; t) 



where the elements of the diffusion matrix are given by 



D 



/ D X x 




■( 




Dpp J 





—ai/mctR \a\ 2 /Ama.R \ 
\a\ 2 /ima R D ) 



(13) 



(14) 



To find a formal solution of (|lj), we use the Fourier representation f(T;t) = J f(T; t) exp [i(xp — px)] dF with f = 
(x,p) T . Eq. ( |T3| ) then leads to 



^fiM = -£ % /(f ;t )-l|D| 
at m 2 



r J D -i r 



/(f;i) 



(15) 



where IDI denotes the determinant of D. The solution takes the form 



/(r ; t) = exp 



--f T G(*)f 

2 v ; 



f(x-pt/m,p;0) 



By substitution into (|15| ) one obtains explicitly the matrix of time-dependent coefficients: 



G(t) = 



D 



D xp + D pp t/2m 



—D 



D pp t/2m D xx - D xp t/m + D pp t 2 /3r 



(16) 



(17) 



Eq. ([l3]) can be interpreted as a Fokker-Planck equation provided the diffusion matrix D is non-negative. Then the 
weight function f(x,p;t) of the pointer states P(x,p) will drift according to the free-particle dynamics. At the same 
time the state of the system will diffuse over the pointer states P(x,p). We can now implement the predictability 
sieve and minimize the production rate for linear entropy by minimizing the width of the Gaussian pointer states. 
In other words, we maximize an under the condition that the diffusion matrix be non-negative. The condition that 
D has a non- negative determinant leads to the condition a|j -I- 2a 2 R a 2 + WDmaj^ai + aj < which can only be 
fulfilled if ai < 0, since an > for (|Io|) to be normalizable. Introducing dimensionless polar coordinates R, <f> by 
a = %/Z?mi?exp(i</>), this condition reads R 2 + 8sin2</> < . The maximum for aji = V DmR cos <ft is reached if the 
equality sign holds, since one could otherwise increase R by holding <f> fixed and thus increase olr. Maximizing an 
under the condition R 2 — —8 sin 2<f> then yields for a the following value a s distinguished by the predictability sieve: 



,1/4 



(V3- 



(18) 



which coincides, up to a small deviation in the numerical coefficients, with the value (||) following from the criterion 
of Hilbert-Schmidt robustness. This above slight departure of a s might be related to the fact that the given form 
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of predictability sieve predicts a degenerate diffusion matrix D. We think, however, that the emerging incoherent 
uncertainties due to the pointer state diffusion must be made proportional to the quantum uncertainties already 
present in the pointer states itself. The "robustness of uncertainties" demands that the matrix C ( pi] ) of quantum 
correlations be proportional to the diffusion matrix D (jf4|) of the corresponding classical coordinates for the pointer. 
From the condition that C = const x D we then obtain again the standard value (^) for a, while C = m/2D x D. 

We shall now discuss our last method to determine the pointer basis, which will involve quantum state diffusion. As 
we see from ( |i2| ) and (|i3|), the quantum state of the system, when expanded as a mixture of pointer states, performes 
diffusion after the decoherence time has elapsed. This diffusion will, by construction, preserve the shape (^) of the 
Gaussian wave packet, and only its center will walk randomly, ft is then natural to ask, whether there is a generic 
QSD process which, first, applies to generic pure initial states and, second, tends to the above specific diffusion process 
for t > t D . 

As is well known, the Fokker-Planck equation ( |l3| ) is equivalent to the Ito-Langevin equation 01| 

dr = Vdt + dX , (f 9) 

where V = (p/m, 0), and dX = (d£, d7r) is the increment of a zero-mean Gaussian white noise with correlation matrix 
Ddt. In case of phase-space diffusion the use of the Ito-Langevin formalism instead of the Fokker-Planck formalism 
is a matter of taste. But the diffusion of the corresponding pointer states ?/>r would be quite awkward in the Fokker- 
Planck formalism. We thus choose the Ito-formalism and apply to the Gaussian pointer states (|Io|). This leads 
to, substituting (x) = x and (p) = p, 

dt{j = --^—fipdt - — {x - (x)f ipdt + {x- (x))ipdz , (20) 

where the index T has been skipped. The deterministic part of the evolution is governed, up to normalization, by 
the same nonlinear wave equation (|?]) which we had obtained from the Hilbert-Schmidt robustness, while the random 
part is driven by the complex Gaussian white noise 

ct 

dz = -d£ + idTr . (21) 



Since the correlation matrix of dX = (d£, d7r) is Ddt, ( |l4] , pi[ ) yields 

M[dzdz*\ = Ddt (22) 

for the mean of the Hermitian correlation, independent of a. But the correlation M [dzdz] still depends on a. A 
most remarkable feature of (^0|) is that any reference to the phase-space variables V — (x,p) has been cancelled. For 
this reason we have omitted the subscript T from ip and extend the validity of the equation to arbitrary initial state 
vectors. It is possible to prove |f6|, [f9[ [2(], |2f| that, starting from whatever initial state i/j(0), the random solution 
tp(t) will tend to be the Gaussian pointer state V'r(t)) where T(t) is governed by the diffusion process fll9|). Eq. fl2C| ) 
is called the Ito-Schrodinger equation of QSD. 

Since the free parameter a still appears in M [dzdz] , we are left with the non-uniqueness problem of the QSD 
equations. If one, however, chooses the distinguished value (g) of a, one finds, using (|2l|), the simple result 

M [dzdz] = M[dz*dz*] = , (23) 

distinguishing a unique QSD. Historically, this unique QSD was in the Fokker-Planck formalism singled out by certain 
invariance considerations |^2|, ^3|. The Ito-Schrodinger equation (20) with the complex Gaussian white noise (22 2^) 



has become the dominating formalism of standard QSD theory |7j extended for arbitrary Markovian reduced dynamics. 
Applying exact forms of robustness criteria we have thus obtained a unique QSD which leads to stationary Gaussian 
pointer states for t ^> to, whose centers undergo a diffusion process. With heuristic forms of robustness one could 
have chosen other QSD equations like in |f9| (an = 2^/Dm) or [|l6) (an = VDm/2\/2). The recent proposal of 
"maximal survival probability" from |24j differs from our first method and does not lead to ( p3| ) [ p"5| . 

In conclusion, we have demonstrated that three different methods of dynamical robustness lead to an essentially 
unique local pointer basis in case of Markovian local dynamics. The corresponding pointer states follow the classical 
trajectories up to a tiny random diffusion. Well-defined stochastic differential equations, known from the theory of 
quantum state diffusion, govern both the formation and the diffusion of pointer states. These states can thus be 
used to characterize local quasiclassical properties. The pointer states are not an absolute property of the system 
in itself, but only characterize certain stability properties with respect to interactions with the environment: They 
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are least sensitive to quantum entanglement, which is why interference terms between them cannot be noticed by 
local observers. They possess thus meaning with respect to an observer-related branch of the total wave function or 
a component corresponding to a potential fundamental collapse |25| , while the interaction with the environment is 
encoded in the choice of our master equation (^). 
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